Global stability criterion for a quantum feedback control process on a single qubit and 
exponential stability in case of perfect detection efficiency 
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Quantum feedback control is a technology which can be used to drive a quantum system into a 
predetermined eigenstate. In this article, sufficient conditions for the experiment parameters of a 
quantum feedback control process of a homodyne QND measurement are given to guarantee feedback 
control of a spin- 1/2 quantum system in case of imperfect detection efficiency. For the case of pure 
states and perfect detection efficiency, time scales of feedback control processes are calculated. 
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I. INTRODUCTION 

In classical control theory, feedback control describes 
processes in which a closed-loop controller is used to steer 
the states or outputs of a dynamical system, which in 
turn effect the inputs of the controller into the system. A 
remarkable approach to feedback control of quantum spin 
systems has recently been elaborated in Here QND 
measurements are utilized to let a quantum system col- 
lapse deterministically onto a predetermined eigenstate. 

In the present article, the stability and the time scale of 
quantum feedback control processes are studied. As a re- 
sult (Theorem[2} , sufficient limits for the experiment con- 
trol parameters are derived to guarantee asymptotically 
stable quantum feedback control processes on a spin-| 
quantum system. It is proved by applying Lyapunov's 
method to the stochastic differential equation govern- 
ing the quantum state evolution, and thus differs from 
the similar result in [l[ proposing numerical methods of 
scmialgcbraic geometry and aiming at applicability for 
higher spin systems where efficient search for Lyapunov 
functions is practically impossible. For the special case 
of pure states and perfect detection efficiency, the quan- 
tum feedback control process is proved to terminate even 
exponentially fast in time. 

The article is organized as follows. First, the notions of 
QND measurements and quantum feedback control for a 
spin-i system are shortly reviewed, before the stochastic 
stability of quantum feedback control processes with im- 
perfect and perfect detection efficiency are studied, and 
the results are shortly discussed. 



stroys, the measured quantum system itself 0, The- 
oretical as well as experimental investigation have been 
intensively made dealing with processes where quantum 
measurements are utilized constructively, for instance 
theoretical considerations of measurement determination 
by the quantum register Q, quantum feedback control 
by continuous measurements [H, 0, 0, Q , especially in 
quantum optics [l|, 0, [ToL [Till , stabilization and purifica- 
tion of two-level systems jl2l . [ill ] , conditional measure- 
ments of coupled quantum dots by a point contact detec- 
tor md or by a SET [MIlEIll, and the conditional 
measurement approach due to Sherman and Kurizki [l9j 
to prepare predetermined field states of atoms trapped 
in optical QED cavities 0, [2l|, HH , as well as a similar 
approach analyzed for spin squeezing in Cs clocks p3| . 

Although these approaches differ considerably in de- 
tail, most of them utilize repeated quantum nondemoli- 
tion (QND) measurements [2J, §3.3], i.e., measurements 
of an observable Y satisfying the self-nondemolition con- 
dition [Y(t),Y(t')] = for all times t, t', as well as 
the back action evasion condition [Y, flint] = 0, where 
flint = J2j \ ® Bj denotes the interaction Hamil- 
tonian between the considered system (the projections 
and the measuring apparatus (Bj). The QND ob- 
servable Y may correspond, for instance, to a Hermitian 
Lindblad operator L, or to a conserved quantity, such as 
a constant of motion of the considered system like polar- 
ization or momentum. 



III. QUANTUM FEEDBACK CONTROL 



II. QND MEASUREMENTS 

In contrast to a measurement in classical physics, a 
quantum measurement inevitably changes, or even de- 
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Due to ideas of Belavkin [f| 0, 0] as well as Wiseman 
and coworkers [1, [25| , an approach to quantum feedback 
control of spin systems has been recently developed by 
van Handel, Stockton, and Mabuchi [l[ . In this approach 
repeated quantum nondemolition measurements are en- 
gineered to let quantum spin systems collapse determin- 
istically onto a previously chosen eigenstate. This idea, 
surprising from a traditional physics perspective, bases 
on the fact that realistic measurements are not instanta- 
neous but take some finite time. If these reduction time 
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scales are of an order attainable by modern digital elec- 
tronics, a quantum filter and a controller can respond 
on the spin system state, feeding back the intermediate 
nondemolition measurement results to a Hamiltonian pa- 
rameter. In this way it is possible, for instance, to deter- 
ministically prepare highly entangled Dicke states [ll[ , to 
generate and utilize squeeze d q uantum states of trapped 
atoms in an optical cavity [231 ] . or to improve quantum 
error correction [26| . 

The quantum stochastic control formalism of van Han- 
del, Stockton, and Mabuchi can be considered as an 
extension of probability theory, and the traditional for- 
mulation of quantum mechanics can be directly recovered 
from it. In [l|, §IV.C] a stabilizing controller is given for 
a quantum system of spin j — ^, schematically depicted 
in Figure [TJ The conditional evolution of the density op- 
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FIG. 1: (Color online) Schema of a quantum feedback control 
process. The QND measurement output y t from the quantum 
system is used to propagate the conditional state of the filter, 
via the feedback signal H(t). The dashed line indicates (clas- 
sical) digital processing, the filter is determined by Eq. (Tljl. 

erator p describing the quantum system depends on the 
probe parameter measurement rate M > in Hz, and 
the detection efficiency r\ £ [0, 1], a pure number. More 
precisely, the conditional evolution of p is determined by 
the stochastic master equation 

dpt = 9*[H(t) t L]ptdt + y/TjJtr[L]ptdWt, (1) 

where H(t) is the control Hamiltonian (with H(t) = 
in case of no feedback), L is an observable one of whose 
cigenstates is the desired final state of the system, Sf* = 
c S*\B(t\L\ is the adjoint generator 

Sf> = -i[H(t), pt] + L Pt L* - \{L*L Pt + p t L*L), (2) 

Jif is the supcropcrator 

Jf[L]p t =L/h + Pt L* - Tr[p t (L + L*)] p t , (3) 

and dWt is the innovations process 

dW t = 2^/M^y t dt-^jTr[p t (L + L*)] dt (4) 

depending on the QND measurement record yt of the out- 
put corresponding to the observable Y(t) [Tl|, §11]. Here 
Y is normalized (i.e., dY" t 2 = dt) and related to L and the 
standard noises A and A* by the Hudson-Parthasarathy 
equation. The innovations dWt is a Wiener increment 
and dW t / dt is a Gaussian white noise. Note that dW t 
is one-dimensional, whereas Eq. ([1]) is operator- valued. 

Usually, the controller Hamiltonian H(t) is determined 
by a few control parameters. It is most desirable if they 
could be adjusted in a way such that the master equation 
has an asymptotically stable fixed point. 



A. Density operator space of a spin-i system 

In the prototypical physical model of homodyne mea- 
surement of a spin system [l[ , where yt denotes the homo- 
dyne measurement record of the output, the observable L 
and the controller Hamiltonian H (t) in Eq. {TJ) are given 
by 



L = VMJ Z and H(t) = B(t)J y 



(5) 



with the usual angular momentum observables J y , J Zl 
and dWt / dt can be identified with the shot-noise of the 
homodyne local oscillator. Here M > is the strength 
of the interaction between the light and the atoms and is 
regulated experimentally by the optical cavity, and the 
control input B{t) is the applied magnetic field. The time 
scale then only depends on the sensitivity 1 / (2y/Mrj) of 
the photodetection per VHz, and the feedback gain pa- 
rameters of the controller B(t). Hence the time scale 
only depends on experimentally controlled parameters. 
For perfect detection efficiency, r\ = 1, the stochastic 
master equation governing a pure quantum system under 
feedback control is one-dimensional and will be tackled 
analytically below. 

For the special case of a quantum system of spin i, 
i.e., a qubit, the space of all density operators p of the 
system is represented by the two-dimensional disc 

D 2 = {(A, i/) 6 R 2 : A 2 + 1/(1/- 1) ^ 0} (6) 

with center (0, \) and radius |, P, §IV.C]. Here the den- 
sity matrix entries are given by pn = v, p 2 2 = 1 ~ v, 
P21 = p*2 = A, and the state (A, v) = (0, 0) to be stabi- 
lized corresponds to p = |1)(1| = diag (0, 1). Note that 
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FIG. 2: (Color online) The density operator space D 2 of a 
spin-i quantum system. The origin = (0,0) corre- 

sponds to the quantum state the point (0, 1) to |0)(0|. 

in this case the imaginary parts of the off-diagonal en- 
tries of p decouple and may be neglected. The stochastic 
master equation ([T]) describing the conditional evolution 
of a single qubit thus is reducible to the two-dimensional 
Ito equation [l[ 



Mr)\t(l-2vt)dW t , 
(7) 



d\t = [B(t)(v t -±)-f\ t ]dt 
dv t = -B{t)\ t dt - 2y/Mrjv t (u t - 1) dW t . 

Its infinitesimal generator [2?], §7.3] is given as 

X=[B{\v){v - |) - f A ]^ - B(X, 



2Mr] 



i>\- 



+ v 2 {v - l) 2 



(8) 



dE[/(s t )] 
dt 
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IV. STOCHASTIC STABILITY OF QUANTUM 
FEEDBACK CONTROL PROCESSES 



A. Stochastic stability 

In control theory stabilization of nonlinear systems is 
usually investigated using Lyapunov theory. In the 1960s, 
the stochastic counterpart of Lyapunov theory [28[ was 
developed by Has'minskii and others. To prove the cen- 
tral results of this article, we first have to define asymp- 
totic stability of stochastic processes. 

Definition 1 Let Wt be a Wiener process on the canon- 
ical Wiener space (fi, J^", P), and let x obey the ltd equa- 
tion on R™, 



dx t = b(x t )dt + a(xt)dWt, 



(9) 



where b, a : R™ — > R™ satisfy the usual growth and Lip- 
schitz conditions for existence and uniqueness of solutions 
[27j |. Then an equlibrium solution of Eq. ((9|), i.e., a 
solution satisfying b{x„) = o~(x„) = 0, is called stable in 
probability if 



lim 

Xq — >X* 



sup \x t 

t>0 



xJ > e 







Ve > 0. 



(10) 



It is called asymptotically stable if it is stable in proba- 
bility and 



lim 

x — >x, 



lim \xt 

t — >oo 







1. 



(11) 



It is called globally stable if it is stable in probability and 



lim \xt 

t — >oo 



xJ = 



= 1. 



(12) 



IS called exponentially stable in p-th moment, p G N, 
[29l | if there exists a pair of constants a, a > such that 



Ellxt-x^] ^ aE[\x Q -x*\ p ] 



(13) 



for all t ^ 0. Especially for p = 1, x* is then called 
exponentially stable in mean, and for p = 2 exponentially 
stable in mean square. The smallest possible value of the 
constant a is referred to as the growth constant, and the 
largest possible value of a as the rate constant or rate of 
convergence. □ 

The first two notions are local properties, whereas the 
third one is a global property of the system. 



B. Imperfect detection efficiency 

If quantum feedback control is performed with only 
imperfect detection efficiency, i.e., < ?y < 1, the follow- 
ing theorem yields a sufficient condition for the global 
asymptotical stability of its final state. Although a simi- 
lar result for a special controller has been shown already 
i n [J §IV.F] by numerical semialgebraic methods, here 
a general analytical criterion relating the controller pa- 
rameters is given. 



Theorem 2 Consider a quantum feedback control pro- 
cess of a spin-l} quantum system, described by the den- 
sity operators in the state space D 2 and by the stochastic 
master equation with the probe parameter measure- 
ment rate M > 0, the detection efficiency r\ G (0, 1] and 
the controller 



B(X, v) = gi\ + g 2 v. 



(14) 



Assume that for the feedback gain parameters gi and g 2 
there exist real constants c > 1, < d < 2(c— 1) such that 
the maximum / max of the auxiliary function f : [0, 5] x 

[-7T, 7r] -> R, 

f(r,6) = [M - (c - l).9i](l - cos 0) - c^r sin 0(1 - cos (9) 
+ [dgi (1 + cos 6) r - (c - l)g 2 - d(gi 2 +M) ] sin 
+ 2dg 2 (l + cos6)(2rcos6 - \) 
- 4Mr)(l - cos0)((l + cos 6») r - \ f 
-4Mr/(l + cos0)((l+cos0)r-l) 2 , (15) 

is negative, i.e., / max < 0. Then the state (0,0) G D 2 , 
corresponding to po = |1)(1|> of the quantum system un- 
dergoing a quantum feedback control process is globally 
stable. Moreover, condition h!5\) implies 



(1 - n)M 4M V 
91 > : ^ 0, — 

c — 1 a 

Proof. Defining the function 

V(X, v)=cv + d\v - A 2 



< 92 < 0. 



(16) 



(17) 



we have V{0, 0) = and V(X, v) > for (A, v) G D 2 \{0}. 
Since ^ A 2 ^ u(l — v) and A ^ — |, we have V = 
(c - l)v + d\v - (A 2 + v{v - 1)) ^ (c -l)v - \v, i.e., 

V{\,v)^{c-\-\)v, (18) 



and especially V > for v > 0. Moreover by (|14]) . 

J?V = ( 9l X + 92v)\(d(v 2 - I - A 2 ) - (c - 1)A] 
- 4fA(di/-2A) 

-4Mr? [A 2 (i/-±) 2 + ^>-l) 2 ] 
= [M - <%iA + g 2 v) - (c - l)<7i]A 2 + d</ 2 i/ 2 (i/ - 5) 
+ [4 ffl ( v -|)-( c _l) fl2 _^£]Ai/ 
-4Af?7[A 2 (^-i) 2 + i/ 2 (^-l) 2 ] . (19) 

We see immediately that J£V(0, 0) = 0. We will prove 
next that JSfV(A 3 i>) < for (A, v) G D 2 \ {0}. Using the 
coordinates (r, 0), where r G (0, |] and G (— 7r, 7t), given 



by 



tan I = -, we have (A, 1/) = r(sin0, 1 



cos 0) . Thus we obtain 



f(r,6) 



&V(\(r,0),i>(r,0)) 
r 2 (l + cos0) 



(20) 



By the assumption of the Theorem, f(r, 8) < for 
(r, 0) G [0, ±] x (-7T, tt), hence ^fV(A, v)<0onD 2 \ {0}. 
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Therefore, V is a strict Lyapunov function on D 2 with 
the only asymptotically stable state (A, v) = 0. 

Since f(r,±w) = 2((1 - rj)M - (c - l) 5 i), as well 
as /(r,0) = -32Mr?(r 2 - (1 + 
assumption l|15p implies (JTHJ). 



16 Mr] /' 
□ 



Therefore, a quantum feedback control process sat- 
isfying the assumptions of Theorem [5] drives a spin-^ 
quantum system to the state (°) £ D 2 , no matter in 
which quantum state the system is initially. This holds 
true even for the worst case, when the initial state is 
Pi = |0)(0| = (°) € D 2 . Note that a standard state re- 
duction measurement would leave the quantum system 
in this state with certainty. 

The next example shows that there indeed exist pa- 
rameter constellations satisfying Theorem [2] 



where t denotes the time duration since the process start. 
In particular, £ D 2 is globally exponentially stable in 



Proof. For a quantum feedback control process of a spin- 
i system, performed with perfect efficiency, the space of 
all density operators p of pure states of the system can 
be reduced to the circle 

S 1 = {(sin0,l + cos0) eR 2 :9e (-ir,n]} (23) 

with center (0, |) and radius |, 0, §IV.C], and the 
stochastic master equation (JTJ) describing the conditional 
evolution of a single qubit, with perfect efficiency T] = 1, 
then is reducible to the one-dimensional Ito equation 

&6 t = (B(t) - Msin0 i cos0 t )d*-A/Msin0 t dW'i. (24) 



Example 3 For the experimentally controlled paramc- 

and the constants 



ters 5i = ™, 52 = ~x> and 7 ? 




FIG. 3: (Color online) The auxiliary function h in (|15[) for 

c = 4, d = 2. Left figure: The 
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92 



4 ' 



graph of f(r,9). Right figure: The corresponding 2-manifold 
imbedded in 3-space (A, v, z), parametrized by A = rsind, 
v = (l + cos<9)r, z = i^fi- for < r ^ i, \0\ < it. 

c = 4, d = 2, the function h in (fl~5|) is negative on D 2 \{0}, 
which may be seen graphically (Fig. [3]) . □ 



C. Pure states and perfect detection efficiency 

For pure quantum states and perfect detection effi- 
ciency, we arc able to estimate the expected running time 
for a quantum feedback control process more precisely, as 
is shown in the following theorem. 

Theorem 4 Consider a quantum feedback control pro- 
cess of a spin-^ quantum system of pure quantum states, 
described by the stochastic master equation (0) with the 
probe parameter measurement rate M > 0, the perfect de- 
tection efficiency r] = 1, and the controller B{t) of 
with feedback gain parameters 

9l > 0, g 2 < 0. (21) 

Then this process has an expected running time 

T qic (M, 9l )=0(c-^+ M W 2 ), (22) 



Choosing the controller as 

B(i) = |-sin0 t + ^(l + co S t ), (25) 

the system stabilizes the state 9 = ±7r, which we mark 
as |1)(1|. Eq. ([M]) has a unique solution on the interval 
[— 7r, 7r] because its coefficients satisfy the sufficient Lip- 
schitz and growth conditions [13, Theor. 5.2.1]. Since the 
diffusion coefficient 



a(0) 



Msin6> 



(26) 



vanishes at 8 = 0, we have to consider the two intervals 
,/~ = (— 7r, 0) and J + = (0, tt) separately. Since moreover 
the drift coefficient 



M 

b(6) = B{6) - — sin 6 cos 6 



(27) 



satisfies b(—n) = b(n) = and 6(0) = g 2 < 0, the state 
6 = is a reflecting barrier ('entrance boundary') for 
states in J~ , but an absorbing barrier ('exit boundary') 
for states in J + , whereas 9 = ±7r both are absorbing 
barriers [10, §5.2.1]. 

To estimate the expected time that a given pure state 
requires to reach the desired state |1)(1|, represented by 
8 = ±7r, we have to compute the expected first exit time 
T for the random variable 6 to leave the interval J~ or 
J + , respectively T is the solution of the inhomogeneous 
linear differential equation [10, §5.2.7], [HI, §10.9] 
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3T 1 2/a . d 2 T ,. n . 8T 
on under the boundary conditions 



T(±7r,t) = T(Q+,t) - Qf{0-,t) = 0, 
and T(0,O) = f(ff). With the change of variable 



1 + cos 6 9 

x = : — 7. — = cot - 

sin0 2 



(28) 



(29) 



(30) 
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we obtain the relations sin 9 = j^rz and cos0 - 



• at — 2 A/ 
YX> -> 2 ~ (1+x 2 ) 
2 da 



^ and b=^ + ¥£. 



M(x 



1+x 2 



(1+x 2 ) 2 



x 2 + l ' 

More- 



over, d9 



d_ 



l+x 

1+x 2 d 
2 dx 



7, i.e. 



i) 2 
,1(1- 



(31) 



Ly n and X n y n . Although the above arguments hold true 
only for |x| < 1, A„ is the eigenvalue corresponding to 
y n for the entire domains of definition, x £ (— oo, 0) and 
(0, oo), respectively. Thus the smallest eigenvalue of L is 



Ai = |j + 1. Now, L can be expressed as Ly 



Hpv' 



with p{x) = exp(— J -^fi dx), 



i.e. 



Hence Eq. ([28)1 is rewritten as ^ 
, d 2 



L 



X 'dx^ + h{x) d-x 



fLT = -1, where 
d 



with 



h(x) 



51 
M~ 



92 9 x — ^X 3 
:I 1 



M 



1 



(32) 



(33) 



and where the boundary conditions T(0, t) = 0, 
T(x,t) — > as ie — > oo, and J^T(x,t) — > as x — > — oo 
hold. Since for |x| < 1 we have — = E^C-l) 1 ^ 2 ", 
we can write h as 



h{x) 



M' 



92 2 

— x 

M 



3a: 3 ) V (34) 



for |x| < 1. An eigenvalue A of L is given by the equation 
Ly = Ay. We use the ansatz y{x) = a k x k . By the 
boundary conditions, y(0) = 0, hence a = 0. Moreover, 
V'( x ) = J2T bakX 1 *' 1 , and y"(x) = H k - l)a k x k ~ 2 , 
hence 



Ly(x) = ^2 k ( k ~ l ) a kX k + h(x) ^2 k 



IdkX 



fc-i 



fc=2 



fe=i 



For |x| < 1, the series X)o°(~ ^Y x2v is absolutely con- 
vergent and we have (Y^i ka k x>c ~ 1 )(Yl'o } (~ ^) u x 2v ) = 
J2?c k with c fc = ka k j: k J-l) k -'x 3k - 2v - 1 , i.e., ci = 



oi(l 



c 2 = 2a 2 (x 



=a 1 > 

x J + x°), c 3 



3a 3 (x 2 



x 6 — x 8 ). Hence Ly(x) has the following coefficients for 
the first powers of x, 



p{x) 



l + x 2 

i 1 4i 

x *f 



(39) 



and r(x) = "^2 ; satisfying the boundary conditions 



y(0,t) = 0, y(x,i) 



as x 



30 



' 9; 



as 



x — ► — oo, and y(x, 0) = 4^/(x). Thus, we have a Sturm- 
Liouville problem on each interval (— oo, 0) and (0, oo) 
separately, possessing the eigensolutions y n correspond- 
ing to the eigenvalues A„ of (|38j) , and our initial problem 
has the solutions 



T(0,t): 



G(9,a,t)r{a)f(a) da 



o Jj± 



G{9,a,t-T)r(a)dadT (40) 



with 



G{6, (X,t)=Y J c -(9i+Mn)nt/2 ^ 



\!J 



where 9 = 2 arctan x and a = 2 arctan x' according to 
(f30"|) . and ||?m|| 2 = Jj± ry 2 dx. Hence an arbitrarily given 
initial state 9 , i.e., f(9) = Sg (9), is pushed exponen- 
tially fast in time t to one of the final states 8 = ±tt, or 
x = 0, because 



T(6,t)=r(6 )J2 



y n (0)y n (8 ) fl + Xn _ A 



1 



^2c- Al *r(0 o )^ 



yn\r v A r 

y n (0)y n {0o) 

n=l 



(42) 



2(SL 
\M 

'9i 

.M 



( 
( 
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+ 3 ) a 3 + 2 • —a 2 - 3a! 



(35) 
(36) 
(37) 



if |x| < 1. Especially, the lowest power of h(x) is the 
term (|j + l). Let y n (x) = a k xk f° r n £ N, i.e., 
a k = for k ^ n. Then the lowest power of Ly n (x) is 
the term {% + n) na n x n . Setting Ly n (x) = X n y n (x) and 
comparing the coefficients, we then get the eigenvalue 



i.e., Eq. (f22|) . for the smallest eigenvalue Ai in (|38j) . □ 



Example 5 Consider the the case g\ = M, gi = — 4^- 
given m Q, Fig. 3(b)]. The n the eigenvalues /i of the ex- 
pected time T to set the quantum system into the state 
|1)(1| given by (J3SD are bounded by fx ^ ^Ai = M. 
Therefore, the greater the measurement rate M, the 
greater is the smallest possible eigenvalue of T. □ 



DISCUSSION 



(38) 



corresponding to the function y n . In turn, once the eigen- 
value is specified, the coefficients a„ +1 , a n+ 2, are 
determined recursively by comparing the coefficients of 



In this article, a stability criterion for a quantum feed- 
back control process has been introduced, as well as its 
expected running time in case of perfect detection effi- 
ciency. In Thcorcm[5J a sufficient limit for the experimen- 
tal control parameters leading to globally stable quantum 
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feedback control processes acting on a spin-i quantum 
system are given. The proof consists of the application 
of Lyapunov's method to the stochastic differential equa- 
tion governing the quantum state evolution under feed- 
back control. 

In Theorem |4] it is shown that, for perfect detection 
efficiency, the quantum feedback control process termi- 
nates even exponentially fast in time. The proof bases 
on the power series ansatz y = a^x for the derived 
equation — 4fLT = — 1 determining the expected time 
T, yielding eigenfunctions y n with corresponding positive 
eigenvalues A n by Eq. i|38|). Mathematically, T is the ex- 
pected first exit time. Theorem @] implies that the ex- 
pected running time T q f c of a quantum feedback control 
process does not depend on the probability neither of the 
desired state, nor of the initial state. Former numerical 
investigations indicate that the running time of quantum 
feedback control algorithms, and thus T q f c , is about a 
tenth of the decoherence time [26[ up to the order of the 
decoherence time (32j . 

Thus we are left with the unsatisfactory situation that 



the general case of quantum feedback control with im- 
perfect detection efficiency could not yet be proved to be 
exponentially stable, in contrast to the marginal case of 
pure states and perfect detection efficiency. Of course, 
the fact that a Lyapunov function proving exponential 
stability could not be found does not mean that there 
does not exist any at all. However, by the proof of The- 
orem [5] such a Lyapunov function cannot be of the form 
p7|) . since with x = (*) € D 2 , \x\ 2 = A 2 + v 2 = 2rv = 
2r 2 (l + cos 6), i.e., with Eq. flU}, 

4f> \x\ 2 ^ &V{x) ^ M 2 ^ 0. ( 43 ) 

but V has a linear term such that there does not ex- 
ist a constant a > satisfying J/fV(x) ^ — aV(x). By 
Theorem^ however, a Lyapunov function may exist sat- 
isfying this criterion at least for pure states. Thus for 
future mathematical investigation the question remains: 
Can quantum feedback control, including purification of 
a mixed state, be exponentially stable in general? 
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